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Abstract 

We consider brane cosmology studying the shortest null path on 
the brane for photons, and in the bulk for gravitons. We derive the 
differential equation for the shortest path in the bulk for a 1+4 cos- 
mological metric. The time cost and the redshifts for photons and 
gravitons after traveling their respective path are compared. We con- 
sider some numerical solutions of the shortest path equation, and show 
that there is no shortest path in the bulk for the Randall- Sundrum 
vacuum brane solution, the linear cosmological solution of Binetruy, 
et al for uo = — 1, — |, and for some expanding brane universes. 

PACS numbers:75.20.Hr,71.10.Hf,71.27.+a,71. 55.-1 
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1 Introduction 



The possibility of using extra dimensions in order to explain features related 
to unified field theories has been advocated several decades ago by Kaluza 
and Klein. After a die out for many years such an idea was reestablished 
in the context of supergravity and string theory, especially in the latter, 
where extra dimensions are required in order that the theory is rendered well 
defined. Meanwhile other problems have been posed in the framework of 
unified theories. One of them is the huge hierarchy between the electro-weak 
scale (~ 100 GeV) and the Planck scale (~ 10^^ TeV). One possibility to 
explain that difference is based on the dynamics of supersymmetry, a very 
beautiful idea that has not, unfortunately, rendered due (and ripe) issues. 
In the usual Kaluza-Klein, and also in the modern proposals to deal with 
extra dimensions, while the 1+3 (physical) dimensions open up to infinity, 
the extra dimensions are confined in a region of the size of the Planck length, 
namely ~ 10~^^cm, staying beyond experimental verification, today or in the 
near future. 

However, it has been recently shown that it is possible to explain the hi- 
erarchy between the electro-weak and the Planck scale by dimensional reduc- 
tion without compactifying the extra dimensions. Moreover, the usual 1+3 
dimensional Einstein theory of gravity can be reproduced on the macroscopic 
distance scale This is quite different from the standard approach, in 

which extra dimensions open up at short distances only, whereas above a cer- 
tain length scale, physics is effectively described by 1+3 dimensional theories. 
Our 1+3 dimensional Universe would be a three dimensional brane living in 
a higher dimensional theory, thus displaying a certain number of additional 
dimensions. A further proposal to deal with the additional dimensions is to 
have them compactified in a submilimeter scale, unifying in a natural way 
the electro weak and Planck scales [Q. 

The possibility of relaxing the constraints on the size of the extra di- 
mensions is very appealing. Such is the case of the Randall- Sundrum (RS) 
model ll], 0, where the Universe is 1+4 dimensional and the Standard Model 
fields are localized on a 3-brane embedded in the 4-dimensional space. Only 
gravitational fields can propagate in all four space directions. At the phe- 
nomenological length scale the Kaluza-Klein zero-modes are responsible for 
the well-posed Einstein 1+3 dimensional theory of gravity and the excita- 
tions provide a correction. Due to the "warp factor" of the brane, a mass 
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scale around that of Planck mass corresponds to a TeV mass scale in the 
visible brane. This explains the hierarchy problem. The cosmological conse- 
quence of this model is also under active investigation P-ITH]. The model 



leads to new perspectives in many interesting aspects such as the question 
of the cosmological constant. 

The construction of the brane-universe can be traced to the study of 
Eg X Eg string theory, presumably 11-dimensional, with the field theory 
limit studied in |jl3[, and where matter fields live in 10-dimensional branes 
at the edge of the space-time. The issue of higher dimensionality and its 
consequences for the early universe have been often discussed in the recent 
literature 
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Problems related to higher derivative gravity and on 
the cosmological constant problem have also been studied, besides the 
AdS / CFT correspondence and Cardy formula . 

In spite of the attractive aspects of the model, causality can be violated, as 
first noticed in |]19| and . We have two choices facing this situation. Either 
we accept the viewpoint that true causality should be defined by the null 
geodesies in the 1+4 universe instead of in the 1+3 brane spacetime or we find 
some mechanism to avoid such a violation on the brane. In the first case, the 
violation must be neglectable in low energy experiments, otherwise, it could 
have been already found. The question is whether it could be substantial in 
cosmology. If the answer is positive, it might help solving the well known 
horizon problem as discussed in |jl9[ and In this paper, we consider 



the following problem. Suppose there are two observers A and B on the 
brane. A can send series of photons or gravitons to B in order to establish 
communication (see Fig. [^). According to the brane cosmology, photons 
travel on the brane while gravitons may travel in the bulk. We consider the 
three questions: (i) what is the shortest path for gravitons, and whether it 
is on the brane or in the bulk; (ii) how earlier the gravitons can arrive at 
B; (iii) what is the difference of the redshift for photons and gravitons after 
they arrive at B. 
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Figure 1: Possible two paths for massless signal propagation. Solid curve PQ 
is a null geodesic on the brane S and broken line PR is a null geodesic in the 



bulk (modified from [|T9|) 



2 Preliminaries 



We shall consider a 5-dimensional metric describing brane cosmology. We 
thus set up a 5-dimensional action of the form 

1 
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gR+ / (fxJ-gC 



The constant is related to the Planck mass as Kg = 
dimensional metric is 



PI 



dsl 



(1) 

The 5- 

(2) 



where represents a maximally symmetric 3-metric. The energy-momentum 
appearing in the Einstein equation Gab = 4'^ab is decomposed as 



AB 



Tab + Z 



AB 



(3) 
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where Tab is the energy-momentum tensor of the bulk matter (in the RS 
scenario it comes from the bulk cosmological constant A, that is, Tab = 
—A6^) and Tab corresponds to the matter content on the brane located at 
y = 0. We are interested in the case where the energy- momentum tensor of 
the bulk matter can be expressed as 

Tb = —^diag{-p - a,p - a,p - a,p - a, -a) . (4) 

Here, a is the brane tension in the RS scenario. The energy- density p and the 
pressure p come from the ordinary matter on the brane and are independent 
of the position. Assuming the Z2-symmetry and a = 0, the Einstein equation 
permits the following exact cosmological brane solution (corresponding to 
A = 0, a = 0, 7jfc = 5jk) 

a = ao{t){l + \\y\) 

n = no{t)il + p\y\) (5) 



b = b 



'0 



where bo is constant in time (a redefinition of y renders it to be 1) and no{t) 
is an arbitrary function (a suitable redefinition of t fixes it to be 1). In the 
above. 



A = -^bop (6) 

= Y + ^) boP (7) 

where is related to the 5-dimensional Newton constant by = SvrGs, 
and the matter equation of state is p = up as usually. 
For a; = — 1 we have the inflationary case, 

ao{t) = H = —p = const., (8) 
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while for uj ^ —1, the usual solution arises. 



ao = t^ .^p = ^, ^=^^y (9) 
Remarkably, the exact solution in the RS model can also be obtained 



Note that the parameters pb and pb in |21] are related to the corresponding 
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ones here in this paper by the relations pb = p + cr,pb = p — a. The solution 
can be written in terms of the function 



a{t,y) 



'1 4i^ + pf 



6A 



3C 



3C 



K5((T + p) 2 



6A 

sinh(/i|?/|) 



On - 

1/2 



cosh(yU?/) 



We now construct the remaining function 



(10) 



n{t,y) 



do(t) 



(11) 



As for eq. (33) in we also have 



p+3— (p + p) = 
ao 



Defining 



A 



'A ^ 

6^5 36 



(12) 



(13) 



and assuming A > and p = up, the Friedman equation can be solved in the 
case C = 0, = 0. For A > 0, 



ao{t) = a^pl 



a 



136A2 

For A = 0, which is the case of RS model, 

2. ^9 ' ^ 



cosh(KgA)f:/g) — 1 + ^ sinh(K5At/g) 



(14) 



ao(t) = a^(K5P^) 



(15) 



where a*, p* are constant (the origin of time being chosen so that ao(0) = 0) 
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3 The shortest cut and the redshift 

Equation for the shortest cut. 

We consider the generic metric (0) for 6 = 1. Consider two points, rA 
and on the brane. In general, there are more than one null geodesic 
connecting to in the 1+4 spacetime. The trajectories of photons must 
be on the brane and those of gravitons may be outside as assumed here. We 
consider the shortest path for both photons and gravitons. Since the 3-metric 
is spherically symmetric, we can omit the angular part and just consider the 
problem for 

dsl = -n^{t, y)dt^ + a^{t, y)f{r)dr^ + dy^ (16) 
The photon path is on the brane {n{t,0) = 1), therefore 

- dt^ + al{t)f{r)dr^ = 0, (17) 

which can be immediately integrated as 

mdr' = r . (18) 

The graviton path is defined in terms of the geodesic equation 

-n^{t,y)dt^ + a'^{t,y)f{r)dr^ + dy'^ = . (19) 

We suppose that the path is parameterized hj y = y(t). Thus the relation 
r = r(t) is obtained by 



f{r)dr = / — dt . (20) 

rA JtA a[t,y) 

We are looking for the path for which, ts reaches its minimum when r = rs- 
For this purpose, we consider the general case 



r f^r')dr'= C[y{t),yity,t]dt • (21) 

Jta JtA 



f{r')dr' = / 

/ta JtA 

For an adjacent path y = y(t) + Sy(t), we have 



f{r')dr'= / C[y{t) + 6y{t),y{t) + my,t]dt (22) 



tA 
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therefore we find the usual condition 

ItA 

The problem is transformed into the Euler- Lagrange problem 

ItA 

In our case, 



5tBC[y{tB),y{tB)-M = 5 f"" C[y{t),y{t);t]dt . (23) 

JtA 

rmed into the Euler- Lagrange problem 
6 C[yit),y{ty,t]dt = ■ (24) 

JtA 



C[y(tly{t)-A= ' / , (25) 

a[t,y) 

and we have 

II = -a-\r^-e)-'"y . (26) 
The Euler- Lagrange equation thus reads 

-y + {- + -)y + { y 

an n a aw^ 

a' 

+ {-n^ - nn') = . (27) 

(X 

Prom this equation we can see that the shortest path is on the brane only 
when 

a' 

-n^ - W = , (28) 



9,(^)=0 . (29) 

Further, if there exists a solution, when y reaches its maximum, where 
y = and y < 0, we have 

d 

- y + {-n'^ - nn') ^ . (30) 

(X 

Thus, — nn', i.e. dy{an~^) must be negative at this point. 
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The equation is a very difficult nonlinear ordinary differential equation. 
There is no guarantee for the existence of the required solutions. In order to 
obtain a solution with both two ends on the brane, we can make the Fourier 
expansion 

+°° i-ji 

y{t) = El/'Sin- -(t-tA) , (31) 

1 = 1 - tA 

a{t,y) = A{y) + J2[a'i{y) sin- -(^-^a) 

1=1 ^gB - tA 

+a'i{y) cos— ^{t-tA)] , (32) 

tgB — tA 

n{t,y) = Ar(y) + ^K(y)sin- -(t-t^) 

1=1 ^gB - tA 

+n,%)cos it-tA)] , (33) 

tgB — tA 

and then substitute back into the differential equation to obtain the coef- 
ficients yi- Here tgB is the time when the graviton arrives at r^, which is 
different from the time t^B when the photon arrives at Tb- It should be 
determined self-consistently by the equation 



r/m-=/-' V'M'.^)-^w ,, (34) 

Jta JtA avt,y) 



.B^v?{t,y)-y\t) 

JtA 

once the solution is obtained. 

If we want to find the path for a graviton so that it can reach the farthest 
within a given time interval [^^,^5], we can also use the Euler-Lagrange 
equation. Then the length difference between geodesies for photons and 
gravitons within a given time interval can be evaluated 



f l(r'),r- = r &^hffLit (35) 

Jta JtA a[t,y) 

f/M*'^/"-^ (36) 

JrA JtA ao{tf) 



Photon and graviton redshift. 
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In general, if A sends out massless signals at x'^ and x'X + dxj^, these 
signals will reach B at x'^ and + dx'^. The relation of x'X, x'X + dx^ and 
can be obtained by solving the geodesic equation. Then the 



redshift of the signal is [22 



QooixB) 9of,ixA)dx'X _ goo{xA) dx\ 
\ gooixA) gou{xB)dx''j^ \ gooixs) dx% 



(37) 



For a static metric such as the Schwarzschild case, it can be shown that 
dx\ = dx%, therefore. 



\ gooixs) 

For the time-dependent RW metric we have 



dec 

d^ 



R{x%) 



(38) 



(39) 



in which case the redshift is given by 



RiA) 
R{x%) 



(40) 



Thus, in the geometric-optics limit, the redshifts in the two cases can be 
systematically discussed. 

Here, we consider that another graviton starts traveling from at a later 
time tA + St A- Its shortest path is in general different from the previous one. 
Let us denote it as = y*{t). Then the time when it arrives at will be a 
later time tgB + 5tgB 



TA JtA+stA ait^y*) 



(41) 



Therefore we have the equality 

^n^{t,y)-y^{t) 



'tA 



a{t,y) 



-dt 



t,B+st,B ^n'^{t,y^) -y^^it) 
tA+stA a{t,y*) 



dt 



(42) 
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For infinitesimal dtA and dtB, we have 



dt 



B 



B 



a{t,y) 

Thus, the graviton redshift is given by 

^gB _ doit A 



dt, 



'^v?{t,y)-y\t) 



a{t,y) 



(43) 



1 - y'itn) 



VgA a'0{tB)\ 1 - y'^itA) 

while for the photon we have 

VgB _ ao{tA) 



(44) 



^gA ao{tB) 



(45) 



4 Examples 



RS vacuum solution^ 0. 
In this case 

niy,t) = aiy,t) = e-''\y\ . (46) 

Eq. ( [27| ) turns out to be 

y + ky' = (47) 

It has two possible solutions, one is y = yA = 0, and the other is y = 
yo + k ln(t — to)- The second solution does not meet our requirement because 
it will not end on the brane. So the shortest path must be on the brane. 



This agrees with the conclusion in . 
The linear cosmological solution. 

We first consider the case = — | so that from (0) fi = 0, a{t,y) = t — y, 
A = — J. The equation is 



Let t — y = u, then 



or 



-it-y)y + y + y'-y' -1 = 
uil + ii^ - 2u'^ = , 

1 d r, 2u 

u= — . 



2^2 - ^3 



u 



(48) 
(49) 
(50) 
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Therefore, 



dii f du , , , 

(51) 



2u — V? J u 

-- cu^ . (52) 



u 2 



2-u 

We can obtain the solution (to and c are two integration constants) 



y = t±^(t-Uf^-^ . (53) 

It is obvious that this path can not end on the brane either. Furthermore, 
we consider the case uj = —1, A = /i = const. ao{t) = e^*. So dy{a/n) = 0. 
Therefore the shortest path is on the brane. 
The general linear cosmological solution [§. 
Consider the case u —1 

ao{t) = A = -| , n = w^ , w = 2 + 3lu (54) 

a{t,y) = f-qt''-^ , n{t,y) = l + '^y (55) 
d{t,y) = qt"-^ - q{q - l)t'i-^y , a'{t,y) = -qt"-^ (56) 
H't^y) = ~q^t~'^y , n{t,y) = qujt~^ (57) 

Letting y = tf{t) in (0), we get a nonhnear differential equation 

-[1 + {2qu; - q)f + {q'u^ - 2q'oo^)f - q'u'fWf + ^ 
+ [q+{2ujq^ -q^ + q-quj)f 
+ {q'^uj - q^uj"^ - 2ujq^ + 2ujq^)p 
+ {2qW-qW)f]{tf + f) 
+ [2qu-q) + q'u'f]{tf + ff 

-[q-qiq-l)f]itf + ff 

+ [{-q-quj) + {q'-Aq^-3qW)f 

+ {3q^uj - 6qW - 3q^uj^)f 

+ (-2gV + 3g^-gV)f 

+ {q'>uj^ - q''uJ^)f] = . (58) 
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Solution z(x)for co=0 

lz(i)=moi 



0.5 



■ t ♦ 



-0.5 



-i.5 r 



-2.5 : 



. v(l)=-10 



1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.1 

x=t/t„ 



1.9 2 



Figure 2: Diagram for y ~ 0.3^p. 



The analysis of such a differential equation is beyond our capability. We 
leave it as it stands and pass to a discussion of some simple cases where 
numerical analysis can be performed. 



The case considered by Binetruy et al. ||2T| is that of a 3-brane universe in 
the 5-dimensional space time with a cosmological constant. For an equation 
of state p = ijjp they found explicit solutions which we use in order to study 
the question of the existence of shortcuts. The solutions are very involved, 
and we first disentangle the equations using a MAPLE program, and further 
on numerically solve the differential equations. We shall consider the matter 
dominated [lo = 0) and radiation dominated (c<j = 1/3) cases. 

The solution of the gravity equations reads |^ 
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Solution z(x)for (£t=0 

lz(l)=01 



» 



♦ ♦ » . . 

T ♦ ; • 



-0.5 



-i.5 



-2.5 - 



1 1.1 1.2 1.3 



1.4 1.5 

x=tlt. 



1.6 1.7 1.8 1.9 



Figure 3: The same diagram as before, with y beginning at the brane. 



where 



ao(t) 



1/9 

(60) 



with and p^, constants. 

In addition, p;, and pb are the matter densities on the brane and on the 
bulk respectively. We have to choose these constants, which we do according 
to the course we are using to discuss the possibilities of shortcuts. We choose 
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Solution z(x)for co=0 

lz(l}=Oj 



14 
13 
12 
11 
10 
9 



>^ 5 



I v(l)=10 

1 v(i)=ia" 



1.1 1.2 1.3 1.4 1.5 1.6 1.7 1. 

x=tlt. 



8 1.9 



Figure 4: Same as before, with positive initial velocity 



the parameters according to the discussion in Binetruy et al. 

Ph = Ph + P , (61) 

where p stands for the ordinary energy density in cosmology given by 

p = p.(aoK)~' , g = 3(l + cu) . (62) 

The intrinsic tension of the brane, pA, has to be identified with Newton's 
constant in order to recover the standard cosmology, that is 



SttG 



6 



(63) 



when p <^ Pa- 

Moreover the 5 dimensional coupling constant /t, the 5-dimensional New- 
ton constant ^(5), and the Planck mass M(5) are related by 

= SttGo) = M^-j* . (64) 
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0.9 



0.8 



0.7 



X 0.6 



0.5 : 



0.4 



0.3 



Solution z(x)for (x>=l/3 

l:(i)=moi 



5 6 
x=M*t 



10 



Figure 5: Diagram for y ~ 0.3£p in the radiation dominated case. Notice the 
plateau followed in the case of lowest initial velocity. 



Furthermore, we follow Randall and Sundrum and relate the bulk energy 
density ps and the cosmological constant density pa by 



PB 



6 



(65) 



At this point all constants are defined in terms of the Planck mass, and 
our discussion of the evolution of gravity signs can be established. 

For the matter dominated case, = 0, we experimented using different 
initial conditions. In general, we prefer to start with ?/ 7^ in order to avoid 
any spurious solution in the differential equation, which is rather singular. 
We thus suppose that y starts at the order of the Planck length. Pictures ^ 
to ^ show some results. We have chosen to plot the adimensional function 
z{x) = pyix), where n corresponds to twice Planck mass units Mp and 
X = t/to, to being the present age of the universe. 
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Solution z(x)for 0=1/3 
[z(l)=0] 



0.5 



• v(l)=10~'" 



0.4 



& 0.3 - 



0.2 : 



0.1 



'-• 



1 2 



5 6 
x=M*t 



10 



Figure 6: Same as before, with vanishing initial position with respect to the 
brane. 

Each graph contains a set of curves corresponding to three typical 
velocities, whose values are shown in the legend of each graph, producing 
similar behaviors. In figures |^ and |^ we use negative initial velocities and, 
independently of the chosen initial point the curve decays and escapes, 
never returning to the same brane. In the case of positive initial velocities, 
picture ^ shows three curves from which we can notice that the greater initial 
velocity is, the further away from the brane the object will travel. 

Summarizing, these graphs show that the gravity wave always "tries to 
follow the brane" , since the y coordinate either drops fast to zero, or drives 
away, which means that the final point reached is far from the original brane. 

We thus conjecture, based on these results, that the shortest path is 
inside the brane, being the one followed by light. However, there is cer- 
tainly room for further paths due to the extremely complicated character 
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Solution z(x)for 0=1/3 

lz(l)=Oj 



0.00 



-0.05 : 



-0.10 



-0.15 



-0.20 : 



-0.25 



-0.30 



5 6 
x=M*t 



v(l)=-10" 
1 v(l)=-10° 



10 



Figure 7: Same as before, with negative initial velocity 



of the differential equation involved in the problem. Moreover, there seems 
to be some attractors in the differential equation, which further complicate 
the matter, rendering a possible solution even more obscure, while opening 
further possibilities of shortcuts, especially in cases where the bulk density 
becomes important. 

Such complications actually do not arise in full in the matter dominated 
case, but can be clearly seen in the radiation dominated era. In these cases, 
solutions are shown in figures |] through 0. Again, we have plotted the 
adimensional function z{x), where x = Mp t in this case. 

Pictures |] and ^ show a plateau behavior for low positive initial velocities; 
however, there is a threshold velocity for which the curve decouples and 
escapes to infinity. Picture shows curves for three negative initial velocities. 
Again, the wave tries to follow the brane from a distance depending on the 
initial velocity value as we had seen in matter dominated case. 
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In the radiation dominated era, u; = 3, attractors are more clearly formed. 
Their meaning is not known and in some cases, where we can avoid dropping 
into them using special initial conditions, it is natural to foresee solutions 
which return to the brane after a roundabout in the bulk, although we have 
to stress that no such solution has been found so far. We leave this more 
difficult numerical problem for a future publication. 
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